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Abstract. Analysis of Baker Nunn camera observations of satellites 1959ai over 1032 days; 
1959?? over 792 days; 1960t 2 over 480 days; 1961& over 150 days; and 1961a$i over 54 days 
yielded results greatly improved over those previously reported. This improvement is primarily 
due to the use of much more data and secondarily to various modifications in the method of 
analysis. As indicated by the discrepancies between results from appreciably different orbits, 
the datum shifts obtained have standard deviations of ±4 to ±23 meters. The tessera! har- 
monic of the gravity field most firmly determined appears to be J 4 i, followed by J 2 2 , Jsi, Jtz, 
J 43 , and J 32 . The principal sources of error suggested are the influence of preassigned variances 
on separation of gravitational coefficients having the same periodic effects on an orbit (e.g., 
J 23 and J 42 ) and the holding fixed with respect to each other stations on the same geodetic 
datum, but they do not seem adequate to explain all systematic discrepancies in the results. 
A comparison of gravitational and geometric geoid heights at station positions gave a mean 
equatorial radius of 6,378,196 ±11 meters. 
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Introduction. This paper describes appreci- 
able improvement over the results obtained by 
Kaula [1963a], due both to additional observa- 
tions and to revisions in the method of analysis. 
The discussion in this paper follows closely that 
of Kaula [1963a] and is limited to changes 
therefrom. 

Observations. The precisely reduced Baker 
Nunn camera observations of 1959«i, 195977, an d 
1960 i 2 from launch until the end of 1961, of 


19618a from launch until the middle of 1961, and 
of 1961aSa in the spring of 1962 were analyzed. 
The observations through mid- 1961 have been 
published in the catalogs compiled by Vets et al. 
[1961-1962]. 

No change was made in the methods of selec- 
tion and conversion of observations. The num- 
ber of observations of each satellite used is given 
in Table 1. 

Geometry . The initial station positions used 


TABLE 1 . Satellite Orbit Specifications 


• 



Satellite 



1959ai 

195977 

1960 t t 

1961 5i 

1961a$i 

Epoch 

1959 Feb. 

1959 Sept. 

1960 Sept. 

1961 Feb. 

1962 Mar. 



28.5 

28.5 

22.0 

20.0 

8.5 


Semimajor axis 

1.304585 

1.334500 

1.250057 

1.252779 

1.568136 


Eccentricity 

0.16582 

0.19008 

0.01146 

0.12135 

0.01197 


Inclination 

0.57381 

0.58212 

0.82434 

0.67835 

1.67316 


Argument of perigee 

3.36062 

3.20403 

2.26377 

2.02733 

4.28853 


Longitude of node 

2.52442 

3.48304 

2.28139 

2.76786 

5.71336 


Mean anomaly 

6.00463 

3.82408 

2.72868 

5.96587 

1.51124 


Perigee motion /day 

+0.09181 

+0.08501 

+0.05186 

+0.08315 

-0.01733 


Node motion/day 

-0.06108 

-0.05712 

-0.05413 

—0.06347 

+0.00367 


Max. A/m, cm 2 /g 

0.21 

0.27 

0.27 

15.9 

0.08 


Min. A/m, cm 2 /g 

0.21 

0.04 

0.08 

15.9 

0.02 


Perigee height, km 

560 

510 

1500 

640 

3500 


Number of days 

1032 

792 

480 

150 

54 


Number of observations 

3513 

3034 

2502 

1395 

552 
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TABLE 2. Datura Shifts (in length urats of 0.378165 meters) 


Datum 

Coordinate 

1959 a ! 

1959 1 ? 

1960 1 , 

1961 *i 

1961 a «, 

Weighted 

Mean 

Americas 

Aui 

- 02.5 

- 02.6 

- 03.7 

- 03.8 

- 06.4 

- 03.8 ± 1.0 


Au 2 

- 04.7 

- 05.2 

- 09.6 

- 11.3 

- 04.2 

- 05.1 =fc 0.8 


A m 3 

- 00.9 

- 00.5 

+ 00.6 

- 01.9 

- 00.3 

- 00.4 ± 0.2 

Europe-Africa 

Aui 

+ 06.5 

+ 07.3 

+ 06.6 

+ 11.6 

+ 04.6 

+ 05.8 ± 0.7 

Siberia-India 

AU2 

- 07.8 

- 07.8 

- 09.3 

- 04.1 

- 10.2 

- 08.9 ± 0.5 


Am 3 

+ 02.0 

+ 01.3 

+ 02.2 

- 01.4 

+ 02.1 

+ 01 . 9 + 0.2 

Australia 

Ami 

- 16.3 

- 19.6 

- 19.6 

- 11.2 

- 26.6 

- 17.3 ± 1.5 


Am 2 

+ 09.6 

+ 06.0 

+ 06.7 

+ 07.5 

+ 03.0 

+ 05.2 ± 1.7 


Am 3 

+ 10.8 

+ 14.6 

+ 10.0 

+ 14.7 

+ 09.4 

+ 10.5 ± 0.4 

Japan-Korea- 

Ami 

- 08.9 

- 11.5 

- 08.3 

- 08.5 

- 06.5 

- 08.9 ± 0.5 

Manchuria 

AU2 

+ 04.1 

+ 05.2 

+ 13.0 

+ 08.7 

+ 09.3 

+ 09.4 ± 0.7 


Am 3 

+ 01.4 

+ 00.1 

+ 01.8 

- 00.1 

+ 04.4 

+ 01.5 ± 0.8 

Argentina 

Ami 

+ 35.6 

+ 37.9 

+ 39.9 

+ 50.7 

+ 34.4 

+ 38.3 ± 1.6 


Au 2 

- 03.7 

+ 03.8 

- 02.4 

- 02.1 

+ 00.0 

- 02.3 ± 0.6 


Am 3 

+ 10.0 

+ 07.5 

+ 09.9 

+ 04.2 

- 06.3 

+ 05.7 ± 3.5 

Hawaii 

Ami 

+ 03.3 

+ 01.3 

- 05.2 

+ 00.4 

- 00.3 

- 04.0 ± 1.6 


Am 2 

+ 06.1 

+ 04.5 

+ 16.0 

+ 01.2 

+ 15.2 

+ 09.2 ± 2.8 


Am 3 

- 45.4 

- 48.6 

— 47.7 

- 67.9 

- 25.0 

- 45.5 ± 3.6 


were the solutions given in Table 1 of Kaula 
[1963a] with corrections for errors in the com- 
puted positions of three stations relative to the 
principal datums provided by I. G. Izsak of the 
Smithsonian Institution Astrophysical Observa- 
tory. Corrections to coordinates u u u 2 , and m 3 
in earth radii are listed in sequence (all values 
times 10 -6 ) : 


San Fernando 

+ 5.6 

-5.0 

-8.5 

Naini Tal 

+2.7 

-5.0 

+4.9 

Curagao 

-1.0 

0.0 

+2.8 

The datum shifts 

listed in Table 2 of 

the pres- 


ent paper apply to the starting coordinates in 
column 4 of Table 1 of Kaula [1963a] with the 
local corrections given above. 

Dynamics. The only change made in the 
dynamical aspects of the treatment was to omit 
entirely the long-period and secular perturba- 
tions that are due to lunisolar attraction, radia- 
tion pressure, and drag by a specified atmos- 
pheric model. For the orbital arc lengths of 10 
to 20 days it was found that these effects were 
adequately absorbed by an arbitrary accelera- 
tion in the mean anomaly. Their inclusion made 
little difference in the solutions obtained for tes- 
seral harmonics or station shifts — if anything, 
they may have distorted the results by shifting 
computed satellite directions farther from those 
observed. 



on page 478 of Kaula [1963a] only three were 
employed: preassigning a covariance matrix V - 
for the starting values of parameters, assigning 
higher weight to the across-track than to the' 
along-track component of an observation, and 
using arbitrary polynomials. 

The covariances and variances preassigned 
were identical with those in Kaula [1963a] with 
the exception of C 2 ,, S 21 , which were held fixed. 

The observational variance employed was 
(0.026 sec) 2 time and (9.2 sec) 2 direction 
for 12-day arcs of 1959®! and 1959 t; ; 
(0.047 sec) 2 time and (13.4 sec) 2 direction for 
20-day arcs of 1960t 2 ; (0.146 sec) 2 time and 
(43.8 sec) 2 direction for 10-day arcs of 19618a; 
and (0.047 sec) 2 time and (13.4 sec) 2 direction • 
for 18-day arcs of 1961a8i. The principal cri- 
terion used in determining the observational 
variances was the x* test; i.e., the quantity m 

s = (f T W~’f - z r M 7 'W' 1 f)/(n - p) (1) 

should average 1 for several orbital arcs, where 
f is the vector of observation equation residuals; 

W is the covariance matrix of observations; z 
is the vector of corrections to parameters; M is 
the matrix parameter coefficients in the obser- 
vation equations; n is the number of observa- 
tions; and p is the number of free parameters. 

In forming the covariance matrix W, observa- 
tions in the same pass were treated as having 
the same timing error. 
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TABLE 3. Gravitational Coefficient Solutions 


Multiply all numbers by a scaling factor of 10~ 6 . 


Coefficient* 

1959ax 

1959*7 

1960 12 

1961 5i 

1961a5i 

Weighted Mean 

ACoo 

4.96 

-8.88 

-0.75 

-18.50 

-9.85 

-2.46 ± 2.36 

A(?20 

-0.06 

-0.06 

-0.05 

-0.29 

0.00 

-0.03 ± 0.02 

C 22 

1.30 

1.36 

1.99 

1.80 

2.52 

1.88 ±0.29 

$22 

-1.74 

-0.76 

-1.63 

-0.32 

-0.89 

-1.38 ±0.17 

Cz 0 

0.97 

0.96 

0.98 

1.01 

0.97 

0.97 ± 0.01 

C31 

1.30 

1.62 

1.53 

-0.96 

1.18 

1.52 ± 0.03 

S n 

0.29 

0.99 

-0.10 

-0.34 

0.46 

0.14 ±0.16 

Cz, 

-0.14 

-0.13 

0.29 

2.35 

-0.84 

-0.02 ±0.26 

$32 

0.49 

0.29 

0.38 

-0.16 

0.98 

0.42 ± 0.06 

Czz 

0.36 

1.11 

0.42 

2.36 

1.70 

0.70 ± 0.26 

Szz 

0.83 

1.11 

0.89 

0.43 

-1.33 

0.76 ± 0.29 

C 40 

0.68 

0.67 

0.61 

-0.35 

0.62 

0.67 ± 0.02 

C41 

-0.38 

-0.38 

-0.33 

-0.48 

-1.00 

-0.33 ± 0.01 

$41 

0.43 

0.53 

0.45 

0.39 

-0.45 

0.37 ±0.15 

C42 

-0.10 

-0.10 

0.02 

0.03 

0.47 

0.01 ± 0.02 

$42 

0.52 

0.68 

0.36 

-0.43 

0.06 

0.35 ±0.15 

C43 

0.18 

0.35 

0.50 

0.44 

0.17 

0.17 ±0.02 

$4. 

0.29 

0.11 

-0.00 

0.16 

0.42 

0.41 ± 0.03 

C44 

0.12 

0.01 

-0.20 

0.20 

-0.24 

-0.01 ± 0.08 

$44 

0.11 

0.22 

0.36 

0.29 

0.32 

0.18 ± 0.05 

CsO 

0.02 

0.03 

0.03 

0.01 

0.02 

0.02 ± 0.01 

C51 

-0.14 

-0.02 

-0.01 

-0.63 

t 

-0.13 ±0.02 

$51 

-0.06 

-0.03 

-0.01 

0.23 

t 

-0.01 ± 0.01 

Ceo 

-0.09 

-0.08 

-0.04 

1.10 

-0.10 

-0.09 ±0.02 

Cei 

-0.01 

-0.03 

0.00 

-0.26 

-0.09 

-0.05 ± 0.03 

$«i 

-0.09 

-0.02 

-0.07 

-0.49 

-0.06 

-0.06 ± 0.01 

Ct2 

-0.04 

0.05 

-0.01 

-0.07 

0.05 

0.01 ± 0.01 

$62 

-0.09 

-0.18 

-0.01 

-0.07 

0.01 

-0.02 ±0.03 

C_zz 

-0.02 

-0.10 

0.15 

-0.02 

t 

0.15 ± 0.01 

$63 

-0.12 

-0.01 

-0.08 

-0.06 

t 

-0.08 ± 0.01 

C_u 

-0.00 

0.06 

0.06 

-0.19 

-0.01 

-0.01 ± 0.01 

$64 

-0.06 

-0.09 

-0.42 

-0.30 

0.03 

-0.03 ±0.07 

C70 

0.12 

0.12 

0.07 

0.09 

0.12 

0.12 ± 0.01 


* C nm and S nm are coefficients of spherical harmonic terms kM/r(a/r) n H nvn such that da — 4 x for 

integration over the sphere^ AC_o 0 and AC 20 are corrections to 0.3986032 X 10 21 (1.0 — 0.00108236P 2 ) cgs. 

t. N° determinations of C 6 1, $51, C 6 3, $«3 were made from 1961a 5i because the partial derivatives of the 
orbit with respect to these coefficients were all smaller than the criterion 0.1 n 1 - 2 [Kaula, 1963a]. 


The arc lengths used were chosen after some 
experimentation as giving a reasonable com- 
promise between magnitude of residuals and 
number of observations. 

The use of arbitrary polynomials was held to 
a minimum; i.e., the only one used was a £ 2 vari- 
ation in the mean anomaly. 

In determining the estimated mean value and 
its standard deviation from several orbital arcs 
of the same satellite, the weighting of a particu- 
lar arc was considered to be proportionate to its 
degrees of freedom. The computer program lim- 
ited to fifteen the number of arcs that could be 
combined at a time. In combining the results of 
several sets of fifteen (or fewer) arcs, the weight 


ascribed to the mean of each set was considered 
to be the inverse of its variance, or standard 
deviation squared. 

In order that the final mean and standard 
deviation reflect as much as possible any syste- 
matic differences which were functions of orbital 
specifications, all sets were combined, with in- 
verse- variance weighting, into four groups: 
1959«! and 1959?;, twelve sets; 1960t 2 , two sets; 
19618!, one set; and 1961aSi, one set. The final 
means and standard deviations given in Tables 
2 and 3 are the result of an inverse-variance 
weighted combination of these four group solu- 
tions. However, for most of the variables, the 
standard deviations from combining the four 
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Fig. 1. Vanguard geoid. Geoid heights, in meters, referred to an ellipsoid of flattening 1/298.24, 
determined from observations of satellites, 1959«i and 1959??. 


groups were smaller than the standard devia- 
tions combining all sixteen sets at once, pri- 
marily because the differences between the 1960^ 
mean and the 1959«i and 1959?; mean were 
smaller than the scatter of 1959aii and 1959?; 
solutions about their own mean. 

To avoid the tendency to prejudge the order 
of magnitude of the solution, which is the main 
defect of the preassigned- variance technique, 
some computer experimentation was tried in 
determining the amplitudes of specified periodic 
variations, in place of harmonic coefficients, in 
holding the reference orbit fixed, and in analyz- 
ing residuals. Applying these methods to one 
satellite at a time did not give as good results as 
the preassigned-variance method, to judge by 
the scatter of solutions. To apply them to data 
from more than one satellite simultaneously re- 
quired considerable program revision which did 
not seem worth while because this method has 
been applied extensively by Izsak [1963]. Other 
changes tried and dropped as unnecessary were 
deleting orbital segments for which observations 
are scanty and holding fixed the station shifts 
obtained from the previous analysis of 1960i 2 
observations. Also dropped was the device of 
weighting observations inversely as their den- 
sity with respect to the phase angle (node- 
GST). 


Results . The analysis described above took . 
much time to apply to the large quantity of 
1959oa and 1959?; data. The attempt to combine 
solutions from different sets of arcs was not 
made until this analysis had been completed. 
Consequently, the good agreement shown by 
Tables 2 and 3 between the results from 1960t a 
on the one hand and from 1959ai and 1959?; on 
the other came as a pleasant surprise. The com- 
bination of results is not as good, of course, as 
is suggested by the formal standard deviations 
given in the tables; in particular, the errors in 
difference of position between stations in North 
and South America — or between stations in 
Europe, Africa, and India — which were held 
fixed with respect to each other, are probably 
several times as great as some of the stated 
uncertainties. The good agreement is even more 
marked for the spatial representations given in 
Figures 1 and 2; e.g., for the seven most extreme 
maximums and minimums in the Vanguard geoid 
of Figure 1, there are maximums and minimums 
in the Echo rocket geoid of Figure 2 agreeing 
within 10° in location and within 11 meters in 
magnitude. The degree of independence in these 
solutions is fairly satisfying. The orbits differ 
by 0.23 in inclination, and 016 in eccentricity, 
the arc lengths used differed in a ratio of 5 to 3, 
and the observational weighting differed in a 
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Fig. 2. Echo rocket geoid. Geoid heights, in meters, referred to an ellipsoid of flattening 
1/298.24, determined from observations of satellite 1960t a . 


« ratio of 3 to 2. It would be very desirable, how- 
ever, to obtain comparable series of observations 
of a satellite of much higher inclination. 

The principal sources of systematic error 
likely to be common to satellites 1959«i, 1959*/, 
and 1960t 2 seem to be (1) that the magnitudes 
of the results will be influenced by the preas- 
signed variances and (2) that the relative posi- 
tions of tracking stations on the same geodetic 
datum may be appreciably in error. 

For a parameter whose effects are fairly dis- 
tinct in periodicities, etc., from those of other 
parameters, it is implausible that its preassigned 
variance could cause a correction that is too 
large or of wrong sign, but it might cause a cor- 
rection that is too small. However, the variance 
actually used in the analyses is not the estimated 
squared magnitude of the correction c r 2 (c), but 
rather iW(c), where N is the number of orbital 
arcs in a set. Since N was always between 10 and 
15, this seems to be no more than a mild re- 
straint preventing occasional ill-conditioned arcs 
from obtaining absurdly large corrections be- 
yond the range of linearity. 

Distortion caused by the preassigned vari- 
ances seems most likely to occur in separating 
gravitational coefficients whose principal effects 
are of the same period; i.e., coefficients J nm and 
J k i such that m = l and n — k is even. The 


most prominent set of such coefficients is 
Ji 2 } and J Q 2 , all of which cause semidaily varia- 
tions of argument 2(Q — 0). A way of remov- 
ing some (but not all) of the influence of the 
preassigned covariances would be to assume 
that what we have determined is not the coeffi- 
cients themselves but the amplitudes of semi- 
daily variations in the orbital elements; e.g., for 
the cos 2(Q — 0) term in the variation of the 
inclination 


A i 


di 

dC 22 


+ 


di 

0C 42 


+ 


di 

dc\ 2 


Co 


( 2 ) 


The semimajor axis and the eccentricity have 
no semidaily variation. If we omit the 1961 
and 196 la Si results and assume that the similar 
1959a! and 1959*7 orbits should be combined, we 
have two sets of eight equations for three un- 
knowns. Using values C 22 = 1.315 X 10 -6 , &» = 
-1.473 X 10- 6 , C* = -0.101 X KV^S* = 0.567 
X 10- a , C« = -0.009 X 10' fl , S Q2 = -0.104 
X 10“° for the combined 1959a! and 1959?; solu- 
tion (corresponding to Figure 1) and using 
values from Table 3 for 1960i 2 we get the com- 
puted amplitudes of periodic perturbations in 
columns 6 and 10 of Table 4. Using these ampli- 
tudes as the observation equation constants and 
solving by the rule of minimizing N(dA2?Z) 2 
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TABLE 4. Semidaily Perturbations of Satellite Orbits 


Satellite 

Element 

(El) 

dEl 

dC 2 2 

dEl 

dCi2 

dEl 

dCw 

10 e X Comp 

AEl c 

dEl 

d*S 2 2 

dEl 

dS\z 

dEl 

d*Sg2 

10 6 X Comp 
A El t 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

1959 ori and 

M 

-2.92 

0.36 

5.74 

-3.93 

2.92 

-0.36 

-5.74 

-3.91 

19597/ 

i 

3.62 

-5.89 

3.83 

5.33 

3.62 

-5.89 

3.83 

-9.12 

combined 

Ci> 

1.68 

12.21 

-24.58 

1.20 

-1.68 

-12.21 

24.58 

-7.09 


O 

-5.54 

4.22 

4.88 

-7.76 

5.54 

-4.22 

-4.88 

-11.10 

1960 12 

M 

-6.04 

0.00 

-0.86 

-12.05 

6.04 

0.00 

0.86 

-9.86 


i 

5.48 

-5.08 

-0.34 

10.81 

5.48 

-5.08 

-0.34 

-10.76 


O) 

-2.59 

20.86 

-2.61 

-4.76 

2.59 

-20.86 

2.61 

-11.70 


a 

-5.07 

-3.26 

7.66 

-10.12 

5.07 

3.26 

-7.66 

-7.07 


yields 

C 22 = 1.85 X 10 -8 S 22 = -1.75 X 1(T 8 

C i2 = 0.05 X 10 -8 S 42 = 0.34 X 10~ 8 

c 62 = 0.10 X 10 -8 §02 = -0.22 X 10~ 8 

All the coefficients are increased over the mean 
in Table 3 except S which hints of ill condition- 
ing. However, it looks as though only Cat and Sw 
might have been significantly reduced by the 
preassigned- variance method. 

The assumption made by Kaula [1963a] that 
the relative positions of tracking stations on the 
same datum should be known through the tri- 
angulation networks with a rms error of ±20 
meters or less was based on standard methods 
of estimating triangulation accuracy [ Bomjord , 
1962, pp. 143-159], as is confirmed by the mis- 
closures of large loops of triangulation: (1) 15 
meters in the 4000-km loop around the western 
Mediterranean [ Whitten , 1952] ; (2) less than 
25 meters in the 10,000-km loop around the 
Caribbean [ Fischer , 1959] ; and (3) within 15 
meters for the 10,000-km loop around the Black 
Sea and Caspian Sea, through Turkestan, and 
connecting in northwest India [ [Fischer , 1961]. 

The connections between the three northern 
hemisphere stations of the EASI system are 
closely associated with loops 1 and 3, and the 
connections between the three northern hemi- 
sphere stations of the Am. system are closely 
associated with loop 2. More in doubt are the 
positions of the stations in the southern hemi- 
sphere in Peru and South Africa, which depend 
on long single arcs of triangulation. A test run 
was therefore made on all the 1960^ data, in 


which the stations Arequipa (in Peru) and Oli- 
fantsfontein (in South Africa) were assumed to 
be on separate datums. The results of this test 
corroborated the assumption as to triangulation 
accuracy; the station in Peru moved 24 meters 
with respect to those in North America, while 
the station in South Africa moved 14 meters . 
with respect to those in Eurasia. The changes 
in the _gravitational coefficients were insignifi- * 
cant — C22, from 1.99 to 2.11 X 10 -8 ; S22, from 
— 1.63 to —1.60 X 10 -8 ; C31, from 1.53 to 1.49 X 
10' 8 ; C 41 , from —0.33 to —0.28 x 10“ 8 ; etc. — 
and the maximum effect on any geoid height in 
Figure 2 was 4 meters. 

There still exists the possibility of errors in 
the local connection of tracking stations to the 
triangulation systems, a matter in which better 
standardization of procedures is needed [Kaula, 
19636]. To check this type of error for stations 
on the major datums we calculate the geometric 
geoid heights corresponding to the final posi- 
tions in rectangular coordinates and then com- 
pare these heights with the gravitational geoid 
heights in Figure 3. To estimate the size of dis- 
crepancies to be expected, we have the geoid 
height variance of 1076 m 2 from autocovariance 
analysis of gravimetry [Kaula, 1959] and a 
mean square height of the satellite geoid of 466 
m a , obtained from the sum of the squares of the 
coefficients in Table 3. If the station positions 
and the equatorial radius were correct, the rms 
expected discrepancy between the geometric and 
gravitational geoid heights due to the inability 
of the satellite orbits to pick up the shorter- 
wave variations would be (1076 — 466) 1/2 = 
±25 meters. 

The results of the comparison are shown in 
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Fig. 3. Combined geoid. Geoid heights, in meters, referred to an ellipsoid of flattening 
1/298.24, determined from observation of satellites 1959«i, 1959?;, 196(h 2 , 196l5i, and 1961«$i. 


Table 5. Applying the mean correction of +31 
meters yields a mean equatorial radius of 6,378,- 
196 ±11 meters and a rms discrepancy of ±38 
meters, which implies a rms radial position 
error of (382 — 252) 1/2 = ± 29 meters. Of the 
stations on the major datums, the 60-meter 
discrepancy for San Fernando causes suspicion 
of local connection error; however, there is also 
a 69-meter discrepancy for Villa Dolores, which 
was free to move to its correct position. 


The agreement of the combined solution in 
Figure 3 with astrogeodetic [ Fischer , 1961] and 
gravimetric [Uotila, 1962] solutions is an im- 
provement over that in Kaula [1963a], particu- 
larly in showing a more pronounced negative 
in the western Atlantic. The discrepancies which 
exist may in part be ascribed to the method of 
analysis of the terrestrial data, since the agree- 
ment is appreciably better with the combination 
of astrogeodetic, gravimetric, and satellite zonal 


TABLE 5. Comparison of Geometric and Gravitational Geoid Heights 


Station 

Datum 

Geometric 

Geoid 

Height, 

m 

Gravitational 
Geoid Height, 
m 

Discrepancy 
for 6,378,165 m 
Radius, m 

Discrepancy 
for 6,378,196 m 
Radius, m 

Organ Pass 

Am. 

-4 

-10 

+6 

-25 

Arequipa 


-16 

-8 

-8 

-39 

Curasao 


-25 

-10 

-15 

-46 

Jupiter 


-18 

-9 

-9 

-40 

Olifantsfontein 

EASI 

+22 

+5 

+ 17 

-14 

San Fernando 


+ 117 

+26 

+91 

+60 

Naini Tal 


-17 

-41 

+24 

-7 

Shiraz 


+14 

-18 

+32 

+ 1 

Woomera 

Au. 

+47 

+24 

+23 

-8 

Tokyo 

JKM 

+54 

+4 

+50 

+19 

Villa Dolores 

Ar. 

+ 104 

+4 

+ 100 

+69 

Maui 

H 

+54 

-12 

+66 

+35 


Note. Geoid heights referred to ellipsoid of equatorial radius 6,378,165 meters, flattening 1/298.24. 
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harmonic data of Kaula [1961], especially for 
western Europe. 

The reference flattening of 1/298.24 is used in 
Figures 1, 2, and 3 to facilitate comparison with 
the results of Kaula [1961, 1963a]. The flatten- 
ing equivalent to the solution obtained for AC j* 
is 1/298.28. The J 2 equivalent is 1082.48 X 1CT 6 . 

In conclusion, it can be said that better ex- 
planations are needed for the systematic dis- 
crepancies indicated by Tables 3 and 5. How- 
ever, considering that the observations used 
herein depended on reflected sunlight ; that they 
were all made more than 3 years before the 
minimum of solar activity; and that the orbital 
■specifications are far from ideal, the prospects 
, are bright for extracting more information on 
* the gravitational field from more recent and an- 
ticipated satellites. It will be of particular inter- 
est to push the analysis to a good determination 
of some sixth- or eighth-degree harmonics to 
see whether or not they corroborate other indi- 
cators of a weak upper mantle. 
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